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ABSTRACT. Bovine coronavirus (BCoV) is a pathogen related to 
enteric and respiratory diseases in cattle worldwide. Enteric (BECoV) 
strains of BCoV are predominant in South America, and genetic 
investigations have been conducted to identify its relationship with 
isolates of respiratory origin (BRCoV). In this study, we used a BRCoV 
strain (BR-UEL11) derived from an outbreak of respiratory disease in 
feedlot cattle in southern Brazil, and compared the partial sequence of 
the polymorphic region of Spike (which was detected and sequenced 
by two distinct reverse transcription-polymerase chain reactions) with 
those of other BCoV strains. The phylogenetic relationship of BR- 
UEL11 with Brazilian BCoV, which is associated with calf diarrhea 
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and winter dysentery (enteric, BECoV; respiratory, BRCoV), and 
classical reference prototypes was analyzed. The analysis showed that 
the BRCoV strains from Brazil clustered with a clade that was distinct 
from most isolates associated with calf diarrhea (BECoV) and ancestral 
prototype strains such as Mebus, Nebraska, and LYVB. Furthermore, 
the BRCoV strains from Brazil clustered with a clade that contained 
recent strains associated with winter dysentery, showing 98-99% 
nucleotide identity with those strains. These results suggested that the 
Brazilian BCoV evolved from being solely enteric to a dual enteric and 
respiratory tropic virus. 

Key words: Bovine coronavirus; Bovine respiratory disease; 

Spike glycoprotein; Dual tropism; Phylogeny 


INTRODUCTION 

Coronaviruses cause respiratory and enteric disease in humans, birds, and other 
mammals (Siddell et al., 1983, Woo et al., 2009). The bovine coronavirus (BCoV) is a 
global pathogen that is recognized as the causative agent of enteric disease; however, the 
contribution of this pathogen towards the development of bovine respiratory disease (BRD) 
was underestimated compared to those of viruses that are primary agents of BRD, such as the 
bovine herpesvirus 1, bovine viral diarrhea virus, bovine parainfluenza-3 virus, and bovine 
respiratory syncytial virus (Fulton et al., 2011). Flowever, recent reports have indicated that 
BCoV might also be a major pathogen of BRD (Lathrop et al., 2000; Boileau and Kapil, 2010; 
Flick et al., 2012). Moreover, isolates of BCoV associated with neonatal and adult diarrhea 
are referred to as BECoV (bovine enteric or enteropathogenic coronaviruses), whereas strains 
isolated from nasal secretions and lung of cattle affected by pneumonia have been classified as 
bovine respiratory coronaviruses, BRCoV (Lin et al., 2000; Boileau and Kapil, 2010). 

BCoV is a ubiquitous betacoronavirus within the subfamily Coronavirinae, family 
Coronaviridae (ICTV: http://www.ictvonline.org/virustaxonomy.asp). The virus is enveloped 
and pleomorphic, with a single-stranded positive-sense RNA genome. It has two major outer 
glycoproteins, namely, Spike (S), which is further divided into domains SI and S2 (Cavanagh, 
1983), and hemagglutinin-esterase (FIE) (Hasoksuz et al., 2002). SI is the receptor-binding 
domain responsible for viral attachment and entry into the host cell, whereas S2 is related 
to membrane fusion. Therefore, S glycoproteins and FIE are major neutralizing antigens of 
BCoV (Saif, 1993; Lin et al., 2000). Mutations in the S protein, mostly in the SI domain, 
are associated with variations in the host range, tissue tropism, pathogenicity, and antigenic 
differences among strains (Saif, 2010). 

Worldwide, BCoV is associated with BRD in cattle of all ages, shipping fever, calf 
diarrhea, winter dysentery (Boileau and Kapil, 2010), and dysentery during the warmer 
seasons (Park et al., 2006; Decaro et al., 2008). Consequently, BCoV can simultaneously 
cause enteric and respiratory infections in the same animal (Park et al., 2007), although genetic 
markers for distinguishing between strains with both types of tropism were not identified 
during molecular analysis of the S glycoprotein (Kanno et al., 2007). Phylogenetic analyses 
based on the S-encoding gene sequence was performed in South Korea (Park et al., 2006), 
Japan (Kanno et al., 2007), and the USA (Fulton et al., 2013), which revealed that ancestral 
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reference strains of BCoV that are exclusively associated with enteric manifestations clustered 
within a distinct phylogenetic group compared to the recent isolates associated with enteric 
and respiratory diseases. 

Since the first description of a Brazilian BCoV associated with calf diarrhea, most 
Brazilian strains of enteric origin have now been described and characterized by molecular 
methods (Jerez et al., 2005; Takiuchi et al., 2008; Stipp et al., 2009; Ribeiro et ah, 2016); in 
addition, assessment of major public databases did not reveal any description of BRCoV in 
other South American countries. This study investigates the recent molecular evolution of 
BCoV by addressing the phylogenetic relationship between the enteric and winter dysentery 
strains based on a polymorphic region of the Brazilian BRCoV. 

MATERIAL AND METHODS 

The BRCoV strain (BR-UEL11) was derived from nasopharyngeal swabs collected 
during an outbreak of BRD in feedlot cattle from the State of Parana, southern Brazil. Nucleic 
acid purification was perfonned as described (Boom et al., 1990), and the presence of BCoV 
was confirmed using a semi-nested polymerase chain reaction (PCR) targeting the N gene 
(Takiuchi et ah, 2006). For genetic analysis, a reverse transcription (RT)-PCR assay using 
the SPK5 pair of primers (Takiuchi et ah, 2008) was used to amplify the sequence encoding 
the nucleotides 1256 to 1904 based on the S gene of the Mebus strain (accession number - 
U00735.2), which includes the polymorphic region (nucleotides 1368 - 1776) (Rekik and 
Dea, 1994). All amplicons were visualized on an ethidium bromide-stained 2% agarose gel. 
The PCR products were purified using the GFXTM PCR DNA and gel band purification kit 
(GE Healthcare, Little Chalfont, UK) and sequenced (ABI3500 Genetic Analyzer, Applied 
Biosystems) using the forward and reverse primers with the BigDye™ Tenninator v3.1 cycle 
sequencing kit (Applied Biosystems, Foster City, CA, USA). 

The phylogenetic analyses compared the BRCoV strain, BR-UEL11 (GenBank 
accession No. KP208779), with 67 Brazilian strains of BECoV and BRCoV, which were 
available in the GenBank. Furthermore, reference strains of respiratory origin, such as OK- 
0514-3 (AF058944) (Chouljenko et ah, 1998), 232NS (DQ320763.1), the winter dysentery 
strain KW10, and the classical calf diarrhea prototype enteric strains Mebus (U00735) (Stair et 
ah, 1972), Nebraska (JQ741969), and the Brazilian LYVB strain (AY606205) (Brandao et ah, 
2006) were included in the study. To include all available data (from the Genbank) pertaining 
to the Brazilian BCoV, phylogenetic analysis was performed by evaluating the partial 
sequence of the polymorphic region (nucleotides 1388-1712). The phylogenetic relationship 
was constructed using the MEGA 7.0 software. The nucleotide sequences were aligned with 
Clustal W, and the phylogenetic analysis was performed using maximum likelihood’s Tamura 
3-parameter model and 1,000 bootstrap replications. 

RESULTS 

The phylogenetic analysis of nucleotide sequences from the polymorphic region 
revealed two groups that form the major clades 1 and 2 (Figure 1). Clade 1 consists of the 
BECoV strains from Brazil (USP- 2 to 14) derived from stool samples of calves affected with 
the winter dysentery WDBR-01 strain, and the ancestral calf diarrhea reference strains LYVB, 
Mebus, and Nebraska. A deletion of 18 nucleotides from positions 1577 to 1594 of the S 
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gene, which is exactly at the polymorphic region (Brandao et al., 2006), was present in all the 
Brazilian strains and the reference LYVB strain within clade 1. 


Sample: 


: respiratory BCoV 
: BR-UEL11 



JQ918875.1 WDBR-R-31 + 
JQ918876.1 WDBR-R-577 + 
JQ918874.1 WDBR-R-34 + 
JQ918873.1 WDBR-R-37 + 
JQ918872.1 WDBR-R-S/N + 
JQ918871.1 WDBR-R-573 + 
JQ918870.1 WDBR-R-715 ♦ 
JQ918869.1 WDBR-R-67* 
JQ918868.1 WDBR-R-766 ♦ 
JQ918867.1 WDBR-R-729 ♦ 
JQ918866.1 WDBR-R-749 ♦ 
JQ918864.1 WDBR-R-58 ♦ 
JF795404.1 WDBR-E-573 
JF795403.1 WDBR-E-S/N 
JF795401.1 WDBR-E-766 
JF795400.1 WDBR-E-749 
JF795399.1 WDBR-E-729 
JF795398.1 WDBR-E-725 
JF795397.1 WDBR-E-715 
JF795396.1 WDBR-E-577 
JF795395.1 WDBR-E-58 
JF795394.1 WDBR-E-54 
JF795392.1 WDBR-E-34 
JQ918877.1 WDBR-R-54^ 


I JF795391.1 WDBR-F-20 
1 JF795402.1 WDBR-B-04 

FJ899737.1 WDBR-96 
JF795405.1 WDBR-9146 


JF795406.1 WDBR-96 
JF795407.1 WDBR-841 
JF795408.1 WDBR-923 
JF795409.1 WDBR-850 
JF795410.1 WDBR-9215 
JF795411.1 WDBR-9224 
JF795412.1 WDBR-868 
JF795413.1 WDBR-587 
JF795414.1 WDBR-934 
JF795415.1 WDBR-887 
JF795416.1 WDBR-299 
JF795417.1 WDBR-933 
JQ965170.1 WDBR-938 


- AY255831 USP1 — 


DQ479422.1 BR-UEL2 
DQ479423.1 BR-UEL3 
DQ479421.1 BR-UEL1 
KP208779 BR-UEL11 (2012) ♦ ^ 

-JF345147.1 WDBR-B3 

JF345145.1 WDBR-B1 
JF345146.1 WDBR-B2 
JF345148.1 WDBR-B4 
JF345150.1 WDBR-B6 
JF345151.1 WDBR-B7 
JF345152.1 WDBR-B8 


DQ320763.1 232NS ♦ 
AF058944.1 OK-0514-3 ♦ 

- AY935646.1 KWD10 

Nebraska (NVSL) 


- U00735.2 Mebus 


I-AY606192.1 USP-02 

"I AY606204.1 USP-14 
-AY606197.1 USP-07 


- AY606199.1 USP-09 


AY606203.1 USP-13 
AY353075.1 WDBR-01 
AY606193.1 USP-03 
AY606194.1 USP-04 
AY606195.1 USP-05 
AY606196.1 USP-06 
AY606198.1 USP-08 
AY606200.1 USP-10 
AY606201.1 USP-11 
AY606202.1 USP-12 
AY606205 LYVB 


Clade 2 


Clade 1 


Figure 1. Phylogenetic relationship, based on the partial sequence of the polymorphic region (nucleotides 1388- 
1712), of selected Brazilian strains of bovine coronavirus with prototype strains of BECoV and BRCoV. The 
nucleotide sequences were identified by GenBank accession number and strain name. BRCoV strains are denoted 
by a black diamond and BR-UEL11 with an arrow. The tree was generated using the maximum likelihood method 
with bootstrap values of 1000 replicates. 


Genetics and Molecular Research 16 (2): gmr 16029580 














































Bovine coronavirus dual tropism 


5 


Most of the Brazilian BRCoV strains clustered in clade 2, including the respiratory 
strain BR-UEL11 identified during the outbreak of BRD (Figure 1), and other respiratory strains 
from Brazil. This clade also contained the respiratory prototypes OK-0514-3 and 232NS, as 
well as the winter dysentery KWD10 strain. Additionally, two subclades, 2A and 2B, were 
fonned by viruses that originated from Brazil. Subclade 2B consisted of the respiratory BR- 
UEL 11 strain, winter dysentery strains of enteric origin, and the BR-UEL 1, 2, and 3 strains 
obtained from diarrheic calves (Takiuchi et al., 2008). The subclade 2A consisted of several 
winter dysentery strains of enteric and respiratory origin and USP-1 from calf stool; however, 
no nucleotide deletion gaps, as observed in the USP strains of clade 1, were identified. 

Calculation of the genetic distance among strains revealed that the BR-UEL 11 strain 
shared 99 to 98% nucleotide identity with similar strains from Brazil within clade 2, but only 
91 to 90% identity with strains from Brazil in clade 1, which harbored the 18 nucleotide gap 
in their sequence. Compared to the prototype strains, BR-UEL 11 showed a nucleotide identity 
of 91% with LYVB, 96% with Mebus, 97% with Nebraska, 98% with KW10, and 99% with 
the OK-0514-3 and 232NS respiratory strains. 

DISCUSSION 

This study demonstrated the phylogenetic relationship among BCoV strains identified 
and analyzed in Brazil predominantly from 2003 to 2010. In addition, a sequence of the SI 
polymorphic region of a BRCoV strain identified in Brazil in 2012 was included in the study. 
Previous studies have suggested that genetic changes in the gene encoding the S protein of field 
strains may interfere with the humoral protection elicited by vaccination because vaccines are 
generally manufactured using the old prototype strains as antigens (Flussain et al., 1991; Lee 
et al., 2010; Kanno et al., 2013). Currently, commercial vaccines against respiratory pathogens 
of cattle containing BCoV-derived antigens are lacking in Brazil. Vaccines containing BCoV 
were developed to prevent calf diarrhea, being formulated only in association with other 
enteric pathogens. However, empirical analysis based on molecular diagnosis (data not shown) 
suggests that respiratory strains of BCoV may be significantly prevalent in Brazilian feedlots. 

Older hypotheses posit that BCoV strains have genetically evolved from the ancestral 
enteric viruses to distinct respiratory and enteric strains, which are distant from the old 
prototypes responsible for calf diarrhea and winter dysentery (Hasoksuz et al., 2002). However, 
based on the analysis of the S gene of South Korean isolates, where strains associated with 
dysentery clustered with a BRCoV isolate and were distant from the classic Mebus prototype, 
it was hypothesized that the same strain could cause both respiratory and enteric disease, 
and that BCoV had evolved to possess dual tropism (enteric and respiratory) (Park et al., 
2006). To investigate this possibility, an experiment was performed with the South Korean 
winter dysentery strain KWD3, which showed that calves developed diarrhea without clinical 
manifestation of respiratory disease, whereas histopathology demonstrated mild to severe 
interstitial pneumonia caused by the KWD3 strain (Park et al., 2007). 

Phylogenetic analysis of the Brazilian strains revealed that the respiratory BR- 
UEL 11 and WDBR-R-(x) strains clustered within a distinct clade that was distant from the 
old enteric BCoV prototypes. Clade 1 contained the old prototypes such as LYVB, Mebus, 
and Nebraska, whereas clade 2 contained more recent strains associated with winter dysentery 
and respiratory tract disease (Figure 1). Similar results were obtained upon analyzing the S 
gene of winter dysentery strains in Japan (Kanno et al., 2007), and with isolates of enteric and 
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respiratory origin from the USA (Fulton et al., 2013). These findings support the theory that 
this phylogenetic distribution resulted from the evolution of BCoV strains to possess dual 
enteric and respiratory tropism. 

In a phylogenetic analysis with strains from the USA (Fulton et al., 2013), one clade 
that contained more recent respiratory and winter dysentery strains further diverged into 
subclades, which varied antigenically. These differences were assessed by vims neutralization, 
using convalescent cattle sera raised against all the American subclades that cross-react with 
each other, and the Nebraska prototype strain from a distinct clade (Fulton et al., 2013). 
Thus, the observed antigenic differences suggest that vaccines manufactured with calf 
diarrhea prototypes, such as the Mebus strain, may have reduced efficiency against certain 
strains circulating in the USA. The same phenomenon might occur with certain strains of 
other countries such as Brazil. Flowever, challenge studies in validated animal models are 
required to ascertain the correlation between strain-specific antigenic distances and efficiency 
of vaccine protection in cattle infected with the current Brazilian field strains of BCoV (Fulton 
et al., 2015). 

In conclusion, we have demonstrated that the respiratory and enteric winter dysentery 
strains of the Brazilian bovine coronavirus included in this analysis are distant from the 
ancestral LYVB, Mebus, Nebraska, and old field strains that clustered with prototypes related 
to calf diarrhea. These findings suggest that the BCoV strains of Brazil probably diverged 
from enteric to dual tropism, and currently the same strain is associated with both enteric and 
respiratory diseases. In addition, since the Spike protein is the major neutralizing antigen of 
coronaviruses, the differences in the S gene sequence between the Brazilian strains and the old 
BCoV prototypes highlight the importance of assessing the efficacy of vaccines raised against 
the current field strains in Brazil. 
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